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ABSTRACT 

We present new results from the only 2D multi-group, multi-angle calculations of core-collapse 
supernova evolution. The first set of results from these calculations was published in Ott et al. 
(2008). We have followed a nonrotating and a rapidly rotating 2O-M model for ~400 ms after 
bounce. We show that the radiation fields vary much less with angle than the matter quantities 
in the region of net neutrino heating. This obtains because most neutrinos are emitted from inner 
radiative regions and because the specific intensity is an integral over sources from many angles at 
depth. The latter effect can only be captured by multi-angle transport. We then compute the phase 
relationship between dipolar oscillations in the shock radius and in matter and radiation quantities 
throughout the postshock region. We demonstrate a connection between variations in neutrino flux 
and the hydrodynamical shock oscillations, and use a variant of the Rayleigh test to estimate the 
detectability of these neutrino fluctuations in IceCube and Super-K. Neglecting flavor oscillations, 
fluctuations in our nonrotating model would be detectable to ~10 kpc in IceCube, and a detailed 
power spectrum could be measured out to ^5 kpc. These distances are considerably lower in our 
rapidly rotating model or with significant flavor oscillations. Finally, we measure the impact of rapid 
rotation on detectable neutrino signals. Our rapidly rotating model has strong, species-dependent 
asymmetries in both its peak neutrino flux and its light curves. The peak flux and decline rate show 
pole-equator ratios of up to ~3 and ~2, respectively. 

Subject headings: hydrodynamics - neutrinos - stars: interiors - supernovae: general 



1. INTRODUCTION 

The basic context of a core-collapse supernova is well- 
established: a collapsing, degenerate core releases an 
enormous amount of energy, ~1% of which couples to and 
unbinds the stellar envelope. Despite decades of research, 
the mechanism of this coupling remains obscure. The col- 
lapsing core rebounds at nuclear densities and launches 
a bounce shock, but detailed simulations in spherical 
symmetry show that the shock wave stalls at 100-200 
km and fails to explode the sta r (IRampp fe; Jankall200Ct 
iLiebendorfer et all l200lL [2001 : iThompson et al.1 12003( 1. 
The best 2D simulations confirm this result; much of the 
shock's energy is lost to nuclear dissociation and escaping 
neutr inos, and the shock stalls (Bcthc 1990; Jank aet al.l 
12001 . 

The neutrino mech anism, first proposed by 
IColgate k, White! ((1966D in its prompt form, posits 
a burst of neutri no emission to energize the shock. In 
its delayed form ([Bethe fc Wilson! H98 5.). an imbalance 
between neutrino absorption and emission behind the 
shock deposits the additional required energy over 
several hundred milliseconds. The rate of energy depo- 
sition in this "gain region" depends on the relationship 
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between the neutrino flux, which declines as ~ r~ , and 
the cooling rate HacT ~ T 6 , which generally falls off 
much more quickly. 

Early 2D simulations successfully exploded su- 
pernova prog enitors by t he d ela yed neutrino 
mecha nism flHerant et al.1 Il994t IBurrows et al.1 
1995t I.Tanka fc Mueller! 119961: iFrver fc Hegerl 120001: 



Frver fc Warren! 120021 12004). though more recent calcu- 
lations with better neutrino transport generally fail to 
obtain explosions. An exception is for ~8-9 M@ stars, 
whose steep density gradients reduce the ram pressure 
of infalling matter. The post-shock pressure needed to 
drive an explosion abruptly drops , and a sub-en e rgetic 
wind-driven supern ova follows ((Kitaura et all [2006; 
IBurrows et all [2007h . There are indications that more 
generic 2D explosions may be obtained with a soft 
nuclear equa tion of state, which results 
neutron star (jBruenn e t al. 2010; M arek 



m a compac t 
Janka 2009). 



However, recent experiments (jShlomo et ah 120061) and 
obser vations of massive neutron star s Jpemorest et al.1 
201Q;) favor a stiff er equation of state. iBruenn et al.l also 
find the shock to stall at a much larger radius than 
other groups, f o r rea sons that are not yet clear. Sec 
iNordhaus et al.l 1)20101 ) and references therein for a more 
thorough discussion of the state of 2D simulations. 

The delayed neutrino mechanism requires ~10% of the 
energy emitted by the core in electron and anti-electron 
neutrinos over the first few hundred milliseconds, a few 
10 51 erg, to be deposited in the post-shock material. 
The details of this energy deposition depend both on 
the neutrino-matter coupling and on the hydrodynam- 
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ics, which determines how long individual particles are 
subjected to net neutrino heating. This latter point has 
been shown to be a strong function of dimension , with 
iMurphv fc Burrows! ([20081 ) and INordhaus et~aT1 <|2010h 
demonstrating that explosions require ~30% less neu- 
trino heating in 2D than in ID, and ^15-25% less in 
3D than in 2D. Two-dimensional core-collapse simula- 
tions display powerful, low-mode oscillations of matter 
behind the stalled shock as a result of the Standi ng Ac- 
cretion Shock Instability fSASI. iFoglizzo fe Taggerll2000l; 
Blond in et al.l 120031: IFoglizzo et all [20071: llwakami et all 
20081 : lYamasaki fc Foglizzol 120081 IScheck et al.l 1200817 7 
Recent 3D calculations confirm the presence of these 
oscillations but find their energy to be spread over 
a wide range of modes forbidden i n axisymmetry 
(iFrver fe Youngl 120071: iFernanded [20101 INordhaus et~aLl 
120101 ) . The coupling of neutrinos to this hydrodynami- 
cally unstable matter is central to the supernova prob- 
lem. 

In this work, we present new, hitherto unpublished fea- 
tures of the only 2D multi-group, multi-angle neutrino 
transport calculations of core-collapse supernova evolu- 
tion ever perform ed. Results from these simulations were 
first published by lOtt et all (|2008l ). These authors stud- 
ied one nonrotating and one rapidly rotating model, and 
discussed the overall behavior of the evolution and dif- 
ferences between models evolved using full multi-angle 
transpo rt and tho se using multi-group flux-limited dif- 
fusion. lOtt et aD also studied the angular character of 
the specific intensity, a quantity that may only be calcu- 
lated using multi-angle transport. We extend these re- 
sults, characterizing the spatial variation and coupling of 
the matter and radiation fields. We also detail the spa- 
tial and temporal variations in estimated neutrino sig- 
na ls observable from Eart h. We obta i n sim ilar results 
to iMarek etUI (|2009l ) and iLund et~aT1 ()2010f ). but with 
more detailed neutrino transport and a different statisti- 
c al techniqu e. 

lOtt et all (|2008l ) demonstrated that multi-angle trans- 
port self-consistently captures the transition from diffu- 
sion to free-streaming, eliminating flux-limiter artifacts 
(though introducing artifacts intrin sic to S n in the op- 
tically thin limit; see iCastor] (|2004l )l They also showed 
that the net neutrino heating rates behind the stalled 
shock are higher in their nonrotating model with multi- 
angle transport than with flux-limited diffusion. Though 
the difference range s from ~5- 10% at early times to ^20- 
30% at later times, lOtt et aD do not find their 2D multi- 
angle model to be significan tly closer to explosion. 

This paper complements lOtt et al.l by examining in de- 
tail the coupling between matter and radiation and the 
relationship between spatial and temporal fluctuations 
in hydrodynamic quantities and neutrino spectral energy 
densities. We find that the magnitude of fluctuations is 
substantially lower in the radiation fields than in the lo- 
cal matter in the critical region behind the stalled shock. 
This results both from the more quiescent hydrodynam- 
ics in the regions where the neutrinos were predominantly 
emitted and from the multi-angle character of the spe- 
cific intensity. Because the local radiation field is an 
integral over contributions from many sources at depth, 
fluctuations tend to be smoothed out. This is an impor- 
tant effect that radial, "ray-by-ray" transport methods 
neglect. 



We also extend lOtt et all (|2008l ) to investigate the re- 
lationship between variations in the hydrodynamics at 
depth and observable variat ions in the neutrin o signal. 
Following the calcluations of Ma rek et al.l (l2009fl and the 
detailed post-processing of ILund et al.l (j2010j ). we then 
estimate the detectability of these rapid neutrino fluc- 
tuations in the water Cherenkov detectors IceCube and 
Super-K. We find that fluctuations characteristic of con- 
vective overturn and shock oscillations could be observed 
in our nonrotating model, but are very weak for our 
rapidly rotating progenitor. Finally, we calculate the ef- 
fects of rotation on the anisotropy of the observed neu- 
trino sig nal We confirm the basic effects of rotation de- 
tailed in Wal der et all (|2005l) . whereby centrifugal sup- 
port suppresses radiation in the equatorial plane. We 
also detail the evolution of the pole-equator anisotropy 
and calculate the probability distribution of observable 
fluxes. 

The paper is organized as follows. In SJ2J we describe 
our numerical techniques and precollapse models. In <|3l 
we present the spatial variation of matter and radiation 
fields, highlighting the spatial uniformity of radiation due 
to the location of the decoupling region and the multi- 
angle character of the specific intensity. In $U we dis- 
cuss temporal oscillations of matter and radiation, and 
present the phasing of dipolar oscillations in these quan- 
tities with the shock position. Having established a con- 
nection between hydrodynamic oscillations and neutrino 
fluctuations, we estimate the detectability of these fluc- 
tuations in Sj5] In we address the impact of rapid 
rotation on the detectable signal, finding an orientation- 
dependent flux asymmetry of up to a factor of ^3 at 
early times and light curves that strongly depend on both 
species and viewing angle. We conclude in <j7l 

2. METHODS 

The numerical methods and initial models us ed in this 
work are described in detail in lOtt et aD (|2008| ) and ref- 
erences therein. In this section, we summarize the im- 
portant points and refer the reader to these references 
for a more thorough discussion. 

2.1. Hydrodynamics 

We use the rad iation-hydrodynam ics c ode VULCAN- 
2D, d escribed by iLivne et all (|2004T ) and iBurrows et all 
(|2007l ). The hydrodynamics module solves the Newto- 
nian Euler equations with artificial viscosity in two steps: 
a lagrangian step followed by a remapping onto an ar- 
bitrary axisymmetric grid. We use a polar grid with 
120 polar angles and 230 logarithmically spaced radial 
zones from 20 to 4000 km. To avoid small cells (and pro- 
hibitive time step restrictions) , we transition to a pseudo- 
Cartesian grid in the inner 20 km. We implement New- 
tonian gravity with a grid-based Poisson solver, and we 
use the fini te-tempe rature nuclear equation of state of 
iShen et all (|1998al )bl). 

2.2. Radiative Transfer 

VU1CAN-2D includes two radiative transfer modules: 
multi-group flux-limited diffusion (MGF1D) and a dis- 
crete ordinates Boltzmann solver (S n ). Flux-limited dif- 
fusion solves the equations of radiative transfer in the dif- 
fusion approximation, using a flux-limiter to handle the 
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transition to free-streaming. The S n solver discretizes 
the polar and azimuthal angles 8 and <fi of the specific in- 
tensity, preserving its multi-angle character. It includes 
emission, absorption, and isotropic scattering with the 
transport cross-section a s , related to the total cross- 
section erf. by (7 s = (1 — (cos0))erfi, but neglects 0(v/c) 
terms, including Doppler shifts and neutrino advection. 
We use a constant spacing in cos 8 and a variable number 
of azimuthal angles (f> to tile the sphere as uniformly as 
possible. With eight polar angles (Sg), we use a total 
of 40 (9,<j>) pairs, while S\2 uses 92 and Sie 162 pairs. 
These angles are the possible directions of the specific 
intensity and, thus, of radiation transport. 

For this study we hav e employed a hyb rid of MGFLD 
and S n , as discussed in lOtt et al.l ((2008) . We collapse 
our progenitor models and follow them after bounce for 
the ensuing 160 ms using MGFLD. We then transition to 
S n by freezing the hydrodynamics long enough to allow 
the radiation to reach equilibrium with the matter. We 
continue to use MGFLD in the inner 20 km, where the 
matter is optically thick to all species and energy groups 
and the diffusion approximation is accurate. We then 
restart the hydrodynamics and follow the evolution of 
the core using this hybrid approach, with MGFLD deep 
in the core and S n at radii greater than 20 km. Because 
the computational cost of S n scales as n 2 , we are only 
able to run the full evolutionary calculations with eight 
polar angles. We have used Sie to compute snapshots of 
the radiation field at 160 ms after bounce. 

We perform all radiation transport with three neutrino 
species and sixteen energy groups. The first energy group 
is at 2.5 MeV, and the rest are logarithmically spaced 
from 5 to 220 MeV. While all neutrinos participate in 
neutral current reactions, only v e and v e are subject 
to charged-current interactions. Because of their simi- 
lar (neutral-current) cross-sections, we group v u , v^, v T , 
and V T into a single "species," which we designate V M ." 

2.3. Progenitor Models 

We begi n our calculations wit h two 20- M Q precollapse 
models of iWooslev et al.l (|2002ft . To study the effects of 
rapid rotation, we apply an angular velocity profile to 
one of our models of the form 

where R is the cylindrical radius, A sets the scale of 
differential rotation, and tto is the central angular fre- 
quency. Because of the computational cost of S n , we are 
unable to run a set of models to explore A-D,q space. In- 
stead, we study the effects of very rapid rotation using 
A = 1000 km and J7o = 7r rads -1 , corresponding to an 
initial central period of 2 s. This central period is faster 
than those currently fayored for evolved m assive stars 
(|Maeder fc Mevnetll2000t iHeger et al.|[2005l ). An impor- 
tant exception is for collapsar models, which require very 
fast pre-c ollapse stellar rotation rates of <5 seconds in 
the core (Woosley 1993). Without angular momentum 
loss, our rotation rate would result in a "millisecond- 
period" protoneutron star, signi ficantly faster than es- 
timated pulsar birth spin rates (lEmmering fc Chevalie"r1 
[19891: IFaucher-Giguere fc KaspiH2006t lOtt et al.l 12006). 
Our two models therefore bracket the spin range of plau- 
sible progenitors. 



3. RESULTS: ANGULAR VARIATIONS IN 
RADIATION AND HYDRODYNAMICS 

There is good evidence that it is easier to explode 
a supernova by the delaye d neutrino mechanism i n 
two dimensions than in one (|Murphv fc Burrowd l2008) . 
The extra dimension opens a rich array of instabili- 
ties, increasing the residence time of particles in the 
gain region where they u n dergo net neutr i no he ating 
([Murphy fc Burrowsl[2008[ ). INordhaus etHI ([201(1 have 
recently extended this result to three dimensions, show- 
ing that a 3D explosion requires ~ 15-25% less neutrino 
luminosity than its 2D analogue. The central importance 
of neutrino heating in multiple dimensions highlights the 
need to accurately capture the coupling between matter 
and radiation in the semi-transparent gain region, over 
which the neutrino optical depths range from a few hun- 
dredths to a few tenths. 

Radiative transfer in the core-collapse supernova con- 
text is fundamentally a seven-dimensional problem, with 
six dimensions of phase space and one of time. In addi- 
tion, there are six neutrino species (three particles and 
their antiparticles) . Because the neutrinos are not in 
thermal equilibrium with their surroundings, simple ap- 
proximations to full Boltzmann transport may miss im- 
portant physical effects. The most obvious example is 
the phenomenon of net neutrino heating, which occurs 
predominantly because the ambient radiation field in the 
region behind the stalled shock is much harder than it 
would be in local thermodynamic equilibrium. In this 
section, we present results from our nonrotating model, 
demonstrating the striking spatial uniformity of radia- 
tion relative to the matter fields throughout the gain re- 
gion. The radiation fields in our rotating model vary 
smoothly from pole to equator. Because the rapid rota- 
tion inhibits convection so strongly, there is little small- 
scale variation in either the radiation or the hydrody- 
namics. 

3.1. Variation of Neutrino and Matter Fields with 

Angle 

The post-shock region of a core-collapse supernova 
displays a range of hydrodynamic phenomena, from 
large-scale shock oscillations and convective overturn to 
smaller scale turbulence. While these phenomena di- 
rectly influence the local properties of matter, their ef- 
fects on neutrinos are more subtle and depend on the 
strength of the coupling between matter and radiation. 
In the optically thick core, the diffusion limit obtains, 
and radiation and matter are in local thermodynamic 
equilibrium. In the free-streaming limit, the neutrino 
properties largely reflect those of matter near the ap- 
propriate neutrinospheres where they are emitted. We 
define the neutrinospheres to be the (energy and species- 
dependent) radii r T where the radial optical depths equal 
2. 
3 ■ 

(*(r,0)p(r,0)) n dr=| (2) 

Because of the rapid transition to free-streaming, the 
above picture suggests two relatively distinct regimes. 
The angular distribution of radiation should resemble 
that of the local matter only where the material is opti- 
cally thick. 
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We present the angular dependence of hydro-dynamic 
and radiation quantities in Fig. Q] at radii of 50 km and 
100 km at three epochs after bounce in our nonrotating 
model. These plots clearly indicate the distinction be- 
tween the two regimes discussed above; only the top-left 
panel shows a region where neutrino spectral energy den- 
sities (indicated at representative energies by dot-dashed 
lines) resemble the angular distribution of density p and 
electron fraction Y c . All quantities arc plotted at fixed 
radius and normalized to their average values over 4-7T 
steradians. The neutrino energies used are the root- 
mean-square values for each species at large radius, 



ations, we introduce the fractional angular variance, de- 
fined for a quantity X as 



_ f fe 2 F v (e)de 
jF„(e)de 



1/2 



(3) 



where F v {e) is the neutrino energy flux spectrum. These 
rms energies increase with time (see lOtt et all I2008L 
Fig. 18). We use the energy groups closest to the rms val- 
ues over the first few hundred milliseconds after bounce: 
16, 21, and 27 MeV for v e , V e , and V M ," respectively. 
Because the couple to matter more weakly than 

electron and anti-electron types, their neutrinospheres lie 
deeper in the protoneutron star core. The V M " spectra 
arc therefore harder, reflecting the higher temperatures 
that prevail where they are emitted. 

As the core evolves during the postbounce phase, varia- 
tions with angle increase in all quantities and at all radii. 
This is due to two effects: 1) growth in the vigor of con- 
vective overturn and large-scale shock oscillations, and 2) 
the contraction of the protoneutron star, so that a fixed 
radius moves outward in Lagrangian coordinates. The 
effects of contraction are most clearly seen at a radius of 
50 km between 160 ms (top-left panel of Fig. [T]) and 250 
ms after bounce (top-center panel). At 160 ms, the base 
of the convective zone lies at about 70 km in radius, and 
the top-left panel displays the properties of a radiative 
zone. By 250 ms after bounce, the base of the convective 
zone has sunk to around 45 km. The region sampled by 
the top-center panel is therefore convectively unstable. 
As a result, fractional fluctuations in p and Y e increase 
from <1% at 50 km and 160 ms after bounce to about 
10% at the same radius at 250 ms. 

While p and Y e vary strongly with angle in the con- 
vective regions, up to a factor of ^2 at late times (see 
bottom- right panel of Fig. [IJ , variations in neutrino en- 
ergy densities remain smaller throughout the gain re- 
gion. Even at 400 ms after bounce at 100 km (lower-right 
panel) , where the variation in density reaches a factor of 
two, fractional variations in neutrino energy densities re- 
main <20%. These neutrinos are predominately emitted 
by matter beneath the base of the convective zone, and 
retain the relatively smooth angular distribution charac- 
teristic of that region. As discussed in the next section, 
their angular distribution also becomes smoother with in- 
creasing radius due to the integrated contributions from 
sources at many angles. This effect can only be properly 
captured by full, multi-angle neutrino transport. 

3.2. Radial Dependence of Neutrino and Matter 
Variation 

Figure [T] indicates a sharp distinction between the an- 
gular dependence of matter in radiative and convective 
zones. To examine the radial dependence of angular vari- 



(* 2 M))n 
<X(r,t)>£ 



(4) 



where (X)n denotes the average value of X over 4ir stera- 
dians. With this definition, a\ = if and only if the 
quantity X is uniform in angle. The fractional angular 
deviation, ctx, is then \J cr'i. 

Figure [2] shows the fractional angular deviation of the 
same quantities plotted in Fig. Q] at 160 ms and 250 ms 
after bounce for our nonrotating model. The hydrody- 
namic variables p and Y c clearly show the base of the con- 
vective zone, which sinks from ^70 km at 160 ms to ^45 
km at 250 ms after bounce. The neutrino spectral en- 
ergy densities are relatively uniform throughout the gain 
region, as suggested by Fig. [1] This is especially striking 
at 160 ms after bounce (left panel), which was computed 
with 16 polar angles for neutrino transport. This snap- 
shot shows increasing angular variation in the neutrino 
fields up until roughly the appropriate neutrinospheres 
(Eq. [2]), and then a transition to a regime where the 
radiation becomes more spatially uniform with radius. 
This is especially apparent in high energy electron neu- 
trinos (like the 27-MeV energy group shown here). Be- 
cause of the ~ e 2 dependence of neutrino cross-sections, 
higher energy neutrinos have neutrinospheres at larger 
radii, where the matter begins to convect. After the ra- 
diation decouples, however, multi-angle effects smooth 
its angular distribution. 

Neutrinos in the model at 250 ms after bounce have the 
same behavior, first showing more variation with angle 
and finally becoming smoother past the neutrinospheres. 
However, because only 8 polar angles were used for neu- 
trino transport, S n artifacts begin to appear at ^100-150 
km in radius. In S n , neutrinos can only be transported 
along the n discrete polar angles defined in the solver. In 
regions of low optical depth, this tends t o confine radi - 
atio n to radial rays, as discussed in, e.g.. iCaston (|2004[ ) 
and lLivne et al.l ([20041 ). 

The radiation field tends to become more uniform with 
angle because of the multi-angle character of the spe- 
cific intensity; its value at a point is an integral over 
contributions from many sources (and therefore many 
angles) at depth. Because an observer at large radius 
can effectively "see" emission from an entire hemisphere, 
variations in the properties of radiation near the neu- 
trinospheres tend to average out. As a result, the illu- 
mination of matter in the gain region is more uniform 
than would be inferred with purely radial transport. We 
quantify this effect in Table [TJ comparing the fractional 
angular deviation for neutrino energy densities at their 
neutrinospheres to that at 150 km in our nonrotating S\§ 
snapshot. Table Q] demonstrates that the fractional an- 
gular deviation is a decreasing function of radius, as dis- 
cussed above, but an increasing function of neutrino en- 
ergy. Because higher energy neutrinos decouple at larger 
radii, they interact more with convecting matter. How- 
ever, multi-angle effects smooth out much of this varia- 
tion. In our S\q snapshot, fractional variations in p and 
Y c are ~10-20% throughout the convective zone, while 
fractional variations in the neutrino spectral energy den- 
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Fig. 1. — Polar plots showing the values of p, Y e , and the spectral energy density of each species near its rms energy (-i/ (ej), Eq.[3j for 
the nonrotating model. E c , E a , and E M denote electron, anti-electron, and mu/tau neutrino types, respectively. The radial coordinate is 
proportional to each quantity's value at fixed physical radius, 50 km in the upper plots and 100 km in the lower plots, normalized to its 
average over 4iz steradians. The matter and radiation fields both show more variation with angle at late times, as the shock oscillations 
grow in amplitude. However, the neutrino energy densities always vary less with angle than the thermodynamic quantities p and Y e . 



sities fall to < 2% by 150 km 
Many groups currently 
dial transport 4 because of 
cost 



use "ray-by-ray" ra- 
transp ort'' because of it s considerably lower 
Buras et all [200l IMarek fc Jankal [2009; 



iBruenn et al 



(Buras et al. 



2010). The "ray-by-ray plus" method 
2006T ) performs accurate lateral transport 



in optically thick regions by keeping terms associated 
with lateral neutrino advection and pressure gradi- 
ents. However, this technique omits the angular flux 
terms, which transport neutrinos relative to the gas, 
and may therefore exaggerate the anisotropy of the 
ne utrino distribution at large radius (c.f. Section 2.3.2 
of iBuras et al.l [2006h . Our results with full multi- angle 

4 first introduced into supernova theory by Burrows ct al. ( 1995) 



transport provide a baseline, which could be used to 
calibrate an otherwise ad-hoc coupling of neighboring 
rays in semitransparent regions. Such an approach 
might provide more accurate transport at a minimal 
additional cost over current "ray-by-ray" methods. 

4. THE RELATIONSHIP OF NEUTRINOS TO 
SHOCK OSCILLATIONS 

The most visually striking feature of two-dimensional 
core-collapse simulations is the large-scale oscillation of 
matter behind the stalled shock. In our nonrotating 
model, the shock position has a strong dipolar compo- 
nent that oscillates at ~50 Hz. In our rapidly rotat- 
ing model, the oscillations are much weaker and emerge 
several hundred milliseconds later. The relationship be- 



G 




Fig. 2. — Normalized angular deviation (y (X 2 )n — {X)q/(X)o in a quantity X, Eq. J4), as a function of radius. Here, u e represents the 

spectral energy density (Ei) of v^, E a ofV e , and E„ of Deviations in p and Y e increase sharply at the base of the convective zone 

and remain high throughout the gain region. Deviations in neutrino energy densities decrease beyond the neutrinospheres (denoted by 
thin vertical lines) due to the integrated contribution of sources from many angles at depth. Throughout the convective region, angular 
deviations of energy density in all but the highest energy v e are much lower than deviations in p and Y e . The neutrino energies shown, 
except 27 MeV for u e , are each species' approximate rms energy (^/ (e 2 ) , Eq. J3j). Increases in the angular deviations of Ej at large radius 
are banding artifacts of the discrete angles used in ijn . 



TABLE 1 
Angular Deviation in u e Spectral 
Energy Density 



e v [MeV] 


r T [km] 


a(r T ) a 


o-(150 km) b 


16 


62 


0.020 


0.0068 


21 


69 


0.029 


0.0096 


2V 


78 


0.049 


0.014 


35 


91 


0.074 


0.022 



a Fractional angular deviation in spectral en- 
ergy density (Eq.[4} at the appropriate neutri- 
nosphere (Eq. [2J. 

b Fractional angular deviation in spectral en- 
ergy density (Eq. |4j at 150 km. 



tween the shock position and oscillations in the post- 
shock m aterial is complex and has been discussed at 
len gth inlFoglizzo fc Taggen pOOOl ) . iBlondin et alJpOOl 
and|Foglizzo et al. ( 2007$. Here, we compute the phasing 
of oscillations in hydrodynamic and radiation quantities 
as a function of radius, tracing the response of matter 
and radiation throughout the post-shock region. We thus 
establish the temporal relationship of the observable neu- 
trino emission to the unobservable hydrodynamics in the 
core. 

We seek to isolate the strongest oscillatory components 
of the shock radius and post-shock hydrodynamic and 
radiation quantities. Because these oscillations are dom- 
inated by low-Z (in particular I = 1) modes in 2D, we use 
the dipole moments of all quantities. We define these to 
be the I = 1 spherical harmonic coefficients in axisym- 



metry (m = 0), 
ai(X(r)) 



V3 
2 



X(r,0)cos0sin0d0 (5) 



in a quantity X. We also define the normalized cross cor- 
relation function between two quantities X and Y with 
temporal offset r as 



rrYV , /dt[X(f)-(X)][Y(t-r)-P0] 



(6) 



Here, (X) denotes the temporal average of X, and ax its 
standard deviation. Thus, C(X, X, 0) is the normalized 
autocorrelation with zero offset and equals unity, while 
a negative value of C indicates an anticorrelation. We 
define the delay between X and Y to be the offset r that 
maximizes their cross-correlation function. We then con- 
vert this delay into a phase difference using a periodicity 
of 19.4 ms. This is the period of the autocorrelation func- 
tion of the shock position, C(R s h, R s h, t), and represents 
an average SASI frequency in our nonrotating model. In 
the following analysis, we use the dipole coefficients a\ 
of radiation and hydrodynamic quantities as the inputs 
to Eq. ©. 

In Fig. [3l we show the phase differences between 
ai(R s h) and the dipole components of hydrodynamic and 
radiation quantities in our nonrotating model. The ver- 
tical axis is the phase in degrees by which ai(i? s h) lags 
a given quantity, while the line thickness corresponds to 
the maximum magnitude of the cross-correlation func- 
tion, with thicker lines indicating stronger correlations. 
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The hydrodynamic quantities show gradual phase shifts 
throughout the post-shock region, while the radiation 
quantities display a constant phase past their decoupling 
radii. These phase shifts may be thought of as sonic de- 
lays, reflecting the fact that the post-shock material is 
not moving as a solid body. 

The phase shifts of individual hydrodynamic quanti- 
ties provide insights into the dynamics of the post-shock 
material. The accretion rate is closely related to the ve- 
locity, and its dipole component is slightly more than 
90° out of phase with that of the shock position. At the 
shock radius itself, the phase difference would be exactly 
90°, as with a simple harmonic oscillator. The phase lag 
varies with depth due to the finite sound speed, reaching 
an offset of 180° near 70 km. This delay is also illus- 
trated by the advected quantity Y c . At large radius, the 
electron fraction peaks on the side opposite to that of the 
shock position due to the expansion of deleptonized ma- 
terial. The peak temperature asymmetry at depth leads 
the maximum extent of the shock by ~90°. At large 
radius, the temperature asymmetry reflects the shock 
asymmetry, as the hot post-shock material expands into 
unshocked infalling matter. We have not shown the to- 
tal heating rate integrated over the gain region, which is 
primarily a function of the volume of the gain region. Its 
dipole moment is almost perfectly in phase with that of 
the shock radius. 

Shock oscillations are hydrodynamic phenomena, and 
the relationship of the shock asymmetries to the radia- 
tion field is weaker than that to hydrodynamic quantities. 
As shown in Fig. [3J however, a clear correlation between 
the shock radius and neutrino asymmetries is present, 
particularly in electron neutrinos. The phasing of this 
relationship is determined near 50 km, and is largely due 
to the minority of electron neutrinos emitted from con- 
vective regions. Anti-electron neutrinos, which decou- 
ple almost entirely at greater depth, display a weaker 
cross-correlation with dipolar asymmetries in the shock 
radius. The asymmetry in both species' asymptotic flux 
is nearly in phase with matter outflows at large radius, 
and leads the asymmetry in the shock position by about 
70° (V e ) and 120° (i/ e ). Though the fractional asymme- 
tries are lower in radiation than in hydrodynamic quan- 
tities, Fig. [3J demonstrates that large-scale oscillations 
are largely responsible for the asymmetries that we do 
observe. A detection of neutrino fluctuations in a real 
supernova would provide strong evidence for shock oscil- 
lations. 

5. DETECTABILITY OF RAPID NEUTRINO 
FLUCTUATIONS 

In the previous section we showed that large-scale os- 
cillations of matter behind the stalled shock are corre- 
lated with oscillations in neutrino flux, particularly in 
v e 's. The detectability of these rapid neutrino fluctua- 
tions depends on the fraction of neutrinos emitted from 
convective regions and on the vigor of large-scale over- 
turn and global oscillations. The fluctuations must be 
sufficiently large to be distinguishable from shot noise in 
a smoothly declining signal. Large-scale hydrodynamic 
oscillations impose a periodicity on these fluctuations, 
opening up a variety of observational tests. In this sec- 
tion, we examine the prospects for detecting rapid fluctu- 
ations in the emission of a real supernova. A robust de- 
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Fig. 3. — Phase (in degrees) by which the dipole component a\ of 
the shock position (Eq.0 lags dipolar radiation and hydrodynamic 
quantities. The line thicknesses are proportional to the values of 
the normalized cross-correlation, with thicker lines indicating closer 
relationships. The phase lag between the shock position and the 
mass accretion rate and Y e gradually decreases, reflecting the fact 
that the post-shock region is not moving as a solid body. Temper- 
ature responds to compression and expansion; its maximum asym- 
metry occurs when the shock position is most symmetric (and its 
velocity is least symmetric). The phasing of neutrino emission is 
determined near a radius of 50 km, and its correlation with the 
shock position is largest in i/ e . Electron neutrinos have the most 
emission in convective regions and the strongest link to turbulent 
hydrodynamics . 

tection of periodicity in the neutrino signal might confirm 
the large-scale oscillations predicted by 2D core-collapse 
simulations. 

5.1. Estimating the Signal 

The neutrino signal detected on Earth is the prod- 
uct of the spectral flux and the detector response func- 
tion. For a nearby supernova, the best counting statistics 
will be provided by IceCube and Super-K, both water 
Cherenkov detectors primarily sensitive to anti-electron 
neutrinos. Each has a response approximately propor- 
tional to the square of the neutrino energy, so that the 

POO 

"signal" cx / f £ (6,R,t)e 2 de, (7) 
Jo 

where the number flux / is a function of viewing an- 
gle 9, supernova distance R, and time. The constant of 
proportionality depends on the fiducial volum e of the de- 
tector, 22.5 k t for Super-K (llkeda et al.l l2007ft and 940 kt 
for IceCube (|Kowarik et all I2009T ) ■ Super- K has almost 
no background, while IceCube will have an estimated 
background rate of approximately 1.34 x 10 3 (ms) _1 
(jKowarik et al.l 120091 ) . Designed for exceptionally en- 
ergetic neutrinos, IceCube would be unable to resolve 
individual neutrino energies or trajectories, but would 
observe a supernova burst as an increase above its nor- 
mal background. Super-K has a negligible background 



and would provide approximate neutrino directions and 
energies, but with far fewer events. These two detec- 
tors also have different temporal resolut ions. IceCube 
bins its data into intervals of 1.6384 ms (Kowarik et al. 
2009), while Super-K resolves the r elative arrival tim e of 
individual events to microseconds flke da et al.ll2007h . 

In Fig. 21 we show estimated signals computed by 
Eq. (J7]) for both our nonrotating and our rapidly rotat- 
ing models. The fractional fluctuations in electron and 
anti-electron estimated signals are ~10% in the nonro- 
tating model and are weak functions of viewing angle. 
The magnitude of these fluctuations is much larger than 
the secular change of the signal over a SASI period. In 
contrast, the estimated signal for our rotating model is 
dominated by smoothly declining accretion luminosity 
and cooling. Rapid rotation stabilizes the post-shock 
region against convection and suppresses oscillations in 
the hydrodynamics and neutrinos. While these "signals" 
do display some periodic modulation, the magnitude of 
rapid fluctuations is comparable to the secular change 
over a period. Orientation effects, which we explore in 
detail in Sj6j dominate in our rapidly rotating model. 

The rapid variations in our neutrino "signals" arise 
from a combination of vigorous convection, modulated 
by the SASI, and neutrino emission from the convert- 
ing region. If the power in shock oscillations is spread 
over a wider range of modes (as has been suggested by 
recent 3D simulations), these oscillations might leave a 
weaker imprint on the neutrinos. Such a possibility needs 
to be addressed quantitatively, ideally with full 3D sim- 
ulations. Sophisticated post-pr ocessing of 3D hydrody - 
namic calculations, like those of lNordhaus et all (2010), 
may also provide insight. 

We also caution that our estimated "signals" neglect 
neutrino flavor oscillations. Because they decouple from 
matter beneath the convective layer, V p " show frac- 
tional fluctuations of < 1% even in our nonrotating 
model. Flavor mixing would make it more difficult to de- 
tect rapid n eutrino fluctuation s by diluting the v e with 
these "i/p." iLund et al.l (|2010( ) found that under an as- 
sumption of complete, energy-independent mixing of all 
species, fractional fluctuations in their estimated "sig- 
nals" declined by ~|. We also show detection results for 
this extreme case in order to bracket the range of physi- 
cally plausible neutrino signatures of flux variations, and 
find similar results. 

5.2. Detecting the Fluctuations 

The estimated neutrino signals shown in Fig. 0] are each 
a superposition of a smoothly declining component and 
a smaller, rapidly varying component correlated with the 
hydrodynamics around the protoneutron star. Here, we 
estimate the maximum distance to our nonrotating su- 
pernova model at which these va riations would be de- 
tectable by IceCube and Super-K. ILund et"aTI (|2010D es- 
timated the detectability of rapid fluctuations in their 
models using high-frequency Fourier components. They 
found that, neglecting flavor oscillations, IceCube would 
be able to detect a few high-frequency components in 
their model at 10 kpc. Here, we use the Rayleigh test 
(|Leahv et al.lll983l) , a different and simpler method com- 
mon in radio and X-ray astronomy, and obtain similar 
results. 

The Rayleigh test is a test for periodicity in a discrete 
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Fig. 4. — Estimated V e signals (Eq. (7)1 and fully mixed signals 
l/3F e + 2/3'Vp" of our nonrotating model and F e "signals" of 
our rapidly rotating model. The fractional fluctuations are as high 
as 10% in the nonrotating model, but <2% in the rotating case. 
To the extent that neutrino flavor mixing does occur, smoothly 
declining l V M " will dilute the rapid fluctuations shown here in V e . 



time series. Given an assumed period, each element of 
the series is assigned a phase <f> and thus a unit vector 
in r-tfi space. The Rayleigh statistic 7Z is the normalized 
magnitude of the vector sum of these elements, 



E 



sm c 



E 



cos ( 



(8) 



over a temporal range consisting of an integer number of 
periods. For a constant signal with only Poisson noise, 7Z 
will be be d rawn from a y 2 d istribution with two degrees 
of freedom (jLeahv et al.lll983|) . 

The Rayleigh test, as described above, takes its null hy- 
pothesis to be a uniform signal. In the case of supernova 
neutrino emission, we wish to extend the null hypothe- 
sis to include a smoothly, secularly changing signal (but 
still without rapid fluctuations). We may accomplish this 
with a slight modification to the Rayleigh test. In addi- 
tion to the phase angle 0, we define a second phase angle 
9 = <f> + tt. We take our temporal range for 9 to be an 
integer number of periods, but with its starting and end- 
ing times offset by one-half period from their values with 
</>. In essence, we use the difference between the first and 
last half-periods to calibrate out any secular change. We 
then define our modified Rayleigh statistic 1Z' to be 



W = 



J(E sin ^ _ E sin6, i 

v i=i j=i 

l / N * N e y 

+ ^ f ^2 cos fa ~ E cos di ) 

\ i = l 7 = 1 ' 



(9) 



Because of the identities sin 9 = — sm(9 + tt) and cos ( 
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— cos(# + 7r), TV reduces to the Rayleigh statistic TZ ex- 
cept for the different ranges over which 6i and <fii are 
defined. In the limit of many detections over many peri- 
ods, N<f, — > Ng and TZ' obeys exactly the same statistics 
as TZ but with a more general null hypothesis. 

One may search for periodicity either at a pre- 
determined frequency or by sweeping through parameter 
space. Here, we perform a coarse sweep, sampling pe- 
riods from 5-40 ms in intervals of 1 ms. We choose a 
threshold value of 13.16 for TZ' , which gives a 5% prob- 
ability of a spurious detection. We then calculate the 
fraction of Monte Carlo realizations of signals computed 
at a given distance that exceed this threshold in at least 
one frequency. These fractions are our estimated proba- 
bilities of detecting high-frequency periodicity. 

Figure [5] shows these detection probabilities as a func- 
tion of supernova distance for both the IceCube and 
Super-K neutrino detectors. We estimate the detectabil- 
ity under both the optimistic assumption of no flavor 
mixing and the pessimistic assumption that the flavors 
blend completely. The true signal will almost certainly 
lie somewhere in between, but will require a detailed cal- 
culation with a given mass hierarchy and mixing angle, 
which is beyond the scope of this paper. Even with the 
pessimistic assumption of complete flavor interchange, 
rapid fluctuations in our nonrotating model would be 
detectable by IceCube if the supernova were to occur 
within ^5 kpc. 
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Fig. 5. — The probability of a 2<r detection of at least one rapidly 
varying component of the neutrino luminosity using the modified 
Rayleigh test (Eq. [ij on Monte Carlo realizations of the estimated 
signals in F e and l/3F e + 2/3'V M " shown in Fig. [4] Depending 
on the extent of neutrino mixing, rapid neutrino fluctuations in 
our nonrotating model may be detectable by IceCube as far away 
as ~10 kpc. High frequency fluctuations would be detectable by 
IceCube only within ~3 kpc even under the optimistic assumption 
that no flavor mixing occurs. 

Figure [S] shows the distance at which any rapidly vary- 
ing component may be detected. To gain physical insight 
into the nature of the shock oscillations, we would need 



a more detailed power spectrum of the high-frequency 
variation. Figure [6] shows the same detection probabil- 
ities as plotted for IceCube in Fig. [5j but as a function 
of assumed period at four distances. A detailed mea- 
surement of the power spectrum of neutrino fluctuations, 
with statistically significant detections at several periods, 
would be possible for a supernova near the Galactic cen- 
ter, given the (optimistic) assumption of no flavor mix- 
ing. 
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Fig. 6. — The probability of detecting high-frequency variation in 
the anti-electron neutrino luminosity with IceCube and the modi- 
fied Rayleigh test (Eq. [9jl at a given periodicity. In our nonrotating 
model, neglecting flavor oscillations, details of the SASI should be 
observable in neutrinos out to ~8 kpc. 

6. THE EFFECTS OF ROTATION ON NEUTRINOS 

Much of the previous discussion has focused on the 
consequences of convection, shock oscillations, and multi- 
angle neutrino transport. Our rapidly rotating model 
is fundamentally different in these respects. Convec- 
tion near the equator is strong ly inhibited by centr ifugal 
support (see the discussion in iFrver fc Heeerl 120001 ) . As 
shown in Figs. 19 and 21 of lOtt et all (]2008[ ), convective 
overturn tends to be confined to the polar regions where 
the centrifugal support is weaker, and strong shock oscil- 
lations do not begin to develop until about 400 ms after 
bounce. The rotating core evolves relatively smoothly, 
and is dominated spatially by quadrupole variations in 
the matter distribution and neutrino flux. The small- 
scale variations presented in fJ3] and shock oscillations 
explored in 21 are both strongly suppressed. 

Rapid rotation creates centrifugal support and an 
oblate distribution of matter around the protoneutron 
star. At small radii, where the matter is optically thick, 
the radiation fields are similarly oblate. The neutrino dis- 
tribution transitio ns to a prolate form beyond the oblate 
neutrinospheres (|Kotake et al.l [2003) . Figure [7] shows 
this transition. We measure the prolateness or oblate- 
ness of the density and neutrino energy densities using 
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the quadrupole coefficient of a spherical harmonic decom- 
position. Normalizing by the average value at (spherical) 
radius r, we define 



0.6 



q 2 (r) _ /dny 2 °(n)x(r) 

oo (r) = §dnY °(n)X(r) 



(10) 



for a given quantity X. 

The general forms of the matter density and neutrino 
fluxes, shown in Fig. are generic in rotating cores. The 
scale and extent of the matter oblateness are set by the 
rotation profile. As discussed in f}2] of this paper and in 
IQtt et al.1 (|2008fl . rapid rotation in our model produces 
a relatively large quadrupole moment in the density and 
in the asymptotic neutrino fluxes. Our model rotates 
more rapidly than any presented i n the previous stud - 
ies of I Janka fc Monchmeven (|1989MWalder et aT1(|2005l ). 
and lKotake et al.1 (120031 ). We now examine the effects of 
such rapid rotation, from the point of view of an observer 
meauring the neutrino flux at a single angle. 

Asymmetries in neutrino flux naturally lead a ran- 
domly oriented observer to infer an incorrect luminosity. 
We explore and quantify this effect as a function of neu- 
trino species and of time in our rapidly rotating model, 
comparing the inferred luminosity, 4irr F v , to the actual 
neutrino luminosity, 



L v = 



f v an. 



(11) 



Here, r represents a radius sufficiently large that all neu- 
trino species and energies are well into the free-streaming 
regime. We use r — 150 km for our calculations in order 
to satisfy this requirement, while minimizing the effects 
of S n artifacts ( MM - 

In ij6.ll we compute the angular distribution of neu- 
trino flux and the probability distribution of inferred lu- 
minosities in our Sie snapshot at 160 ms after bounce. 
In £ 36. 2[ we show how these distributions evolve in time. 
Because of the limited angular resolution of our evolu- 
tionary calculation, S n artifacts in the asymptotic fluxes 
become significant at late times. We remove this effect 
at the cost of angular resolution by integrating over each 
of the eight bands in polar angle. 

6.1. Flux Asymmetries: Snapshots 

In Fig. 13 we present results for our 160 ms Sie rotat- 
ing snapshot. The left panel shows the total neutrino 
energy flux in each species as a function of /i = cos 9. 
Large pole-equator flux asymmetries are present in all 
species, ranging from a factor of 2.1 in 77 e to a factor of 
3.2 for "Vfj," This corresponds to an overestimate of the 
true luminosity by 100% (v e ), 60% (V e ) or 110% f'V) 
at the poles, or an underestimate by 27% (^ e ), 23% (77 e ), 
or 34% ("fp") at the equator. These asymmetries are 
comparable to or slightly la rger than those reported by 
Uanka fcM onchmcvcr ( 1989) (who find as much as a fac- 
tor of 3) and lWalder et al . (2005) (who reach a factor of 
2.5). This is likely due to the extreme rotation of our 
model, w hich is higher than in any of the models pre- 
sented by Wal der et al.l ([20051 ). In contrast, pole-equator 
asymmetries are < 4% in the nonrotating model even 
though it manifests strong convective plumes and eddies 
(cf. ©. 
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Fig. 7. — Normalized quadrupole moments (Eq. 1 10 I t of the neu- 
trino energy densities and matter densities in our rapidly rotating 
model at 160 ms after bounce. E c , E a , and E M denote the total 
energy density in v e ,v e , and respectively. The centrifugally- 

supported matter remains oblate throughout the computational 
domain. The radiation fields follow the matter in optically-thick 
regions, but become prolate as they decouple, with neutrinos es- 
caping more freely near the poles. 

In the right panel of Fig. [SJ we convert the angular dis- 
tribution of neutrino flux into a probability distribution 
by assuming the observer to be randomly oriented with 
respect to the axis of rotation. Such an observer would 
have a 34% chance of being within 20% of the actual lu- 
minosity in u e , a 54% chance in 77 e , and a 28% chance 
in 'Vp." These distributions have significant tails, par- 
ticularly for 'V M "s for which there is a 14% chance of 
observing at least 1.5 times the mean flux and a 2.3% 
chance of observing at least twice the mean. 

All of these calculations neglect neutrino flavor oscil- 
lations. Should the species mix, anti-electron neutrinos 
detectable in IceCube and Super-K could blend with mu 
and tau antincutrinos, increasing the observable asym- 
metry in v e . Flavor oscillations would thus increase the 
uncertainty in inferring the true neutrino luminosity. 

6.2. Light Curve Asymmetries 

A core-collapse supernova's neutrino light curve is a 
probe of the physical processes deep within the col- 
lapsed core. There is a burst of v e emission as 
the shock wave reaches the electron neutrinospheres 
([Burrows fc Ma zurck 1983). The early postbounce phase 
is powered by accretion onto the core, while much of the 
total energy emerges during the t ens of seconds dura- 
tion o f the protoneutron star phase ([Burrows fc Lattimerl 
1986). A rapidly rotating supernova will have a light 
curve strongly dependent on viewing angle, particularly 
during the accretion-powered phase a few hundred mil- 
liseconds after bounce. In this section, we compute the 
light curve shape as a function of viewing angle in both 
our nonrotating and our rapidly rotating models. 

Unfortunately, our evolutionary calculations lack the 
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Fig. 8. — Left panel: sphericized neutrino luminosities (Anr 2 F v ) as a function of viewing angle for the Sig snapshots. The short dashed 
lines indicate the true luminosities (integrated over 47r steradians) for the rapidly rotating model. Fluxes are calculated at 150 km, after 
they have nearly reached their asymptotic values, but before S n banding becomes noticeable. Right panel: probability of inferring less 
than a given fraction of the true neutrino luminosity in our rapidly rotating snapshot. Even at 160 ms after bounce, there is a 14% chance 
of overestimating the luminosity in 'V M " by at least 50% and a 28% chance of underestimating it by at least 30%. This is especially 
important if flavor oscillations mix l V M " with the more uniform (and more easily detectable) v e and 77 e . 



polar angle resolution to accurately calculate the angu- 
lar distribution of neutrino fluxes at large radius. As the 
neutrinospheres sink and optical depths decrease at fixed 
radius, S n banding becomes significant at smaller radii. 
We remove these artifacts at the cost of angular resolu- 
tion, by calculating average fluxes over angular wedges of 
22.5° (an S n band, or one-eighth of a hemisphere). We 
may still measure the evolution of the pole-equator asym- 
metry by choosing one band near the pole and a second 
near the equator. Because the flux is nearl y symmet- 
ric wi th respect to the equator (cf. Fig. 18 of IQtt et al.l 
(2008)), we only show results for polar viewing angles in 
the range 0° < 9 < 90°. We normalize all observed fluxes 
to their values at 160 ms after bounce to show evolution- 
ary differences. These are arguably more significant than 
the total inferred power, as they are insensitive to uncer- 
tainties in distance. 

We show our results in Fig. |H1 The nonrotating model 
(left panel) shows little variation in its light curve with 
angle. In this case, variations with angle are < 10%, com- 
parable to or lower in magnitude than rapid temporal 
fluctuations. However, orientation effects dominate the 
light curves of the rapidly rotating model (right panel). 
The shape of the light curve shows little variation with 
angle in V e , but declines twice as much in thermal 'V M "s 
near the equator as near a pole. The flux in v e , while 
declining by ^35% near the equator, actually increases 
for an observer near a pole. Because water Cherenkov 
detectors like IceCube and Super-K are primarily sensi- 
tive to v e , the degree to which the less uniform V M " mix 
with V e will determine not only the observed luminosity, 
but also the shape of the observed neutrino light curve. 



7. CONCLUSIONS 

We have presented new results from the only 2D core- 
collapse calculations with full multi-angle neutrino trans- 
port. We have computed the spatial distributions of mat- 
ter and radiation fields, finding that the radiation fields 
are more uniform than matter fields throughout the cru- 
cial gain region. This is both a result of the location of 
the decoupling region (primarily beneath the base of the 
convective region) and due to the multi-angle character 
of the specific intensity, so that the radiation at a point 
is an integral over many sources at depth. This latter 
effect can only be properly captured by full multi-angle 
neutrino transport. Our calculations therefore provide 
an important calibration test of methods, like the "ray- 
by-ray" approach, that do not naturally handle lateral 
transport. 

We have measured the relationship of the radiation 
fields to the large-scale, low-order mode shock oscilla- 
tions, the most striking feature of 2D core-collapse sim- 
ulations. This provides an important way to connect 
observable fluctuations in the neutrino luminosity to the 
otherwise unobservable hydrodynamics around the col- 
lapsed core. We have also calculated the phasing of 
the dipole moments of hydrodynamic quantities with the 
shock radius, demonstrating both the temporal structure 
of the oscillations and their va riation with depth. 

Following earlier work by iMarek et al.l (|2009f ) and 
ILund et all (pOlOh . we have estimated the detectabil- 
ity of temporal fluctuations in the neutrino luminosity 
by the current water Cherenkov detectors IceCube and 
Super-K. We use a novel variant of the Rayleigh test 
for periodicity. With our modification, the test takes 
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Fig. 9. — Light curves of the nonrotating (left panel) and rapidly rotating (right panel) models over four ranges of polar viewing 
angles from to 90°, normalized to 160 ms after bounce. The angular ranges integrate over S„ artifacts, and we use line thicknesses 
proportional to the solid angle subtended. Because the flux is nearly symmetric about the equator, we only show viewing angles for the 
northern hemisphere. The nonrotating model has orientation effects of < 10%, comparable in magnitude to temporal fluctuations. The 
rapidly rotating model has a nearly constant asymmetry in V e flux, but its "v^" light curve declines twice as fast near the equator as it 
does near the north pole. The v & flux decreases by ~35% near the equator, while actually increasing near a pole. Because of the differences 
between light curves of the various species, flavor oscillations could have a large effect on the shape of observed light curves. 



its null hypothesis to be any smoothly varying signal, 
not simply a constant signal. This modified Rayleigh 
test applies to a wider range of problems in which a pe- 
riodic, rapidly varying component is superimposed on 
any secularly changing signal. We find that, neglecting 
neutrino flavor oscillations, we expect to measure a rea- 
sonably detailed power spectrum with IceCube for our 
nonrotating model within ^8 kpc. Under these same 
(likely optimistic) assumptions, Super-K would be able 
to detect rapid fluctuations out to a distance of only ^3 
kpc. Given the pessimistic assumption that all neutrino 
species mix uniformly, rapid fluctuations in our nonrotat- 
ing model would be detectable in IceCube out to ^5 kpc, 
but a detailed power spectrum would probably lie beyond 
the capabilities of current detectors. The true signal will 
almost certainly lie somewhere in between these extremes 
of no and total mixing. This raises interesting prospects 
for the next Galactic supernova, if the low-mode shock 
oscillations in our nonrotating model are comparable in 
magnitude to those found in Nature. Additional and 
complementary information could be gained by the ob- 
servation of gravitational waves that carry with them the 
imprint of m ulti-dimensional dynamics in core and post- 
shock region (lOttll2009HMurphv et al.ll200a iMarek et al.l 
[200l lYakunin et al.ll2O10( h 

We also analyzed the angular distribution of radiation 
in our rapidly rotating model and its effects on the lu- 
minosities and light curves measured by randomly ori- 
ented observers . We find pole-equator asym metries at 
least as large as lJanka &: Monchmeverl (I1989D. and larger 
than any of the models studied by iWalder et al.1 (|2005l ) 
(though ours also rotates somewhat faster). The asym- 



metries are strongly dependent on neutrino species, with 
"i/p" showing the greatest asymmetry (a factor of ^3 
at 160 ms after bounce) and v e the least (a factor of 
~2). The light curves of our rapidly rotating model are 
also strong functions of species and of angle, and again, 
the "i/jj" show the greatest asymmetry and v e the least. 
By 550 ms after bounce, the pole-equator asymmetry 
remains a factor of ^2 in T> e , but has grown to ~6 in 
V„." The asymmetry in any observed signal thus de- 
pends strongly on the degree of neutrino flavor mixing. 
Because of the strong dependence of the light curves on 
angle, a randomly oriented observer could face consider- 
able uncertainty in the timescale, not simply the magni- 
tude, of the neutrino emission. 

By treating the neutrino-matter coupling accurately 
in the crucial gain region, our calculations have enabled 
a detailed study of the spatial and temporal character 
of the radiation fields. These will enable the calibra- 
tion of less accurate (and much less expensive) radiative 
transfer techniques. They also demonstrate the strik- 
ingly uniform character of radiation in the gain region. 
Our estimates of the detectable features of our models 
also suggest that the next Galactic supernova may of- 
fer an observational test of the dipolar shock oscillations 
common to 2D core-collapse simulations, and highlight 
the uncertainty associated with neutrino measurements 
of a rapidly rotating core. 
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